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1.  Introduction. 


Recently  a  lot  of  work  has  been  done  in  the  area  of  robustness  of  multivariate 
tests.  Robustness  typically  carries  three  meanings:  null,  nonnull,  and  optimal¬ 
ity.  Roughly  speaking,  a  test  is  defined  to  be  null  robust  if  its  null  distribu¬ 
tion  remains  the  same  for  a  class  of  distributions  including  the  distribution 
under  which  the  test  is  derived.  Similarly  we  define  the  nonnull  robustness  of 
a  test  by  the  invariance  of  its  nonnull  distribution  in  a  class  of  distributions 
including  the  underlying  distribution.  Finally,  a  test  is  called  optimality  robust 
if  an  optimal  property  the  test  enjoys  can  be  extended  to  a  class  of  distributions 
including  the  distribution  under  which  the  optimality  holds.  In  a  series  of 
articles  published  by  the  present  authors,  the  various  relationships  connecting 
these  robustness  concepts  have  been  clarified,  tools  for  establishing  robustness 
of  many  multivariate  tests  have  been  developed  and  extensively  used  (see  e.g. 

Kariya  (1977,  1980,  1981a,  1981b),  Kariya  and  Sinha  (1985),  Sinha  (1984), 

Sinha  and  Giri  (1985),  Sinha  and  Das  (.1985)).  It  turns  out  that  in  all  the  cases 
considered  so  far  optimality  robustness  along  with  null  robustness  hold  for  many 
familiar  tests  for  a  broad  class  of  nonnormal  distributions.  However,  it  should 
be  noted  that  in  these  approaches  nonnormality  at  the  sacrifice  of  independence 
has  been  considered.  Therefore,  in  dealing  with  independent  populations,  this 
approach  is  not  quite  natural. 

In  this  paper  keeping  (between)  independence  of  the  two  populations  we  con¬ 
sider  some  testing  problems  involving  location  and  scale  parameters  under  some 
conditions  on  the  underlying  distributions.  To  be  specific,  in  Section  2  we 
have  considered  the  problem  of  testing  the  equality  of  two  location  parameters 
without  any  scale  parameter  and  shown  that  the  test  based  on  the  difference 
between  the  two  sample  means  is  conditional  UMPI,  conditionally  given  two  ancillary 
statistics.  However,  it  is  not  unconditionally  so  unless  either  the  populations 
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are  normal  or  there  is  only  one  observation  from  each  population.  In  Section  3  we 
consider  the  problem  of  testing  the  equality  of  two  scale  parameters  with  or  with¬ 
out  the  presence  of  location  parameters  and  prove  that  the  standard  normal-theory 
UMPI  F-test  is  optimality  robust. 

In  Sections  4  and  5  analogous  results  have  been  obtained  in  association  with 
standard  tests  considered  in  life-testing  models  (exponential  distribution). 

It  should  be  pointed  out  that  in  these  tests  null  robustness  does  not  hold 
mainly  because  of  independence. 

As  a  technical  tool,  Wijsman's  (1967)  representation  theorem  is  used. 


2 .  Tests  of  Equality  of  Two  Location  Parameters  (Non-Normal  Case) . 

*  * 

Let  X  :  (n^+l)xl  and  Y  :  (n2+l)*l  be  independent  random  vectors  with  pdfs 
(with  respect  to  Lebesgue  measure) 


(2.1)  fQ  (x  )»  q1[(x*-01e1) '  (x^-O^)],  fQ  (y*)  *  q2[  (y*-02e2)  '  (y*-02e)  ] 


*  ni+l  *  n2+l  n  +1  n  +1 

where  x  e  R  1  ,y  e  R  ^  ,  e.  «  (1,...,1)’€  R  ,  e  ■  (1,. . . ,1) ’ e  R  2  * 

n  +1  n2+l 

C'l’P2e  R  ql  anc*  ^2  are  P^s  over  R  and  R  respectively.  Making  a  suitable 

orthogonal  transformation,  without  any  loss  of  generality,  we  can  write 

*  #  n2 
x  =  (x,x|)',  y  »  (y,yp  ’  where  x^  e  R  ,  y^e  R  and  assume  that  the  pdfs  are 

(2.2)  fQ  (x*)  -  q1[(x-ei)2+ x^},  fQ  (y*)  -  q2l ^y_92)2  +  yiyl^  * 

The  problem  of  testing  H^:  0^  ■  vs  H^:  9^>  02  remains  invariant  under  the  group 

G  of  transformations,  G  *  G^*  G 2>  where  G^  is  the  group  of  translations: 
x-*-x+c,  y-*-y+c,  -°° <  c<°°,  and  G 2  is  the  group  of  orthogonal  transformations: 

0(np  *  0(n2)  acting  on  x^  and  y^  as:  r2^1*  ^i€  ®^ni^ »  1  ■  1»2. 


Clearly  T  *  (x-y,  x’x  ,  y’y  )  is  a  maximal  invariant.  Using  Wijsman’s  (1967) 
representation  theorem,  the  ratio  of  the  pdfs  of  T  under  nonnull  to  null  is 
obtained  as 


(2.3)  R  *  dP I  /dPT 


j°°  q1[(xfc-91)2  +x^x1]q2[(y+c-02)2  +7^1^ 


V  H0 


,oo 

q1[(x+c-61)2  +x^x1]q2[(y+c-ei)2  +y^y1]dc 
—00  **"  ** 


oo 

f  qifc2+x[*i 


3q2  [  (c-z+d)  +y’y1]dc 


,00 

q1[c2+x’x1]q2[(c-z)2  H-y^ldc 

—00 


z  »  x-y,  6  =  0-6, 


£ 


q  x ( ( c-6 ) 2+*[x1 ] q  2 [ ( c-z ) 2+y  Jy 1 ] dc 


r°° 

I  q1[c2+x’x1]q2[(c-z)2+y^y1]dc 
•  —00  ** 


To  derive  an  optimum  invariate  test  based  on  R,  note  that,  conditionally  given 
the  two  ancillary  statistics  xjx^  and  yjy^,  the  ratio  of  the  conditional  pdfs 
of  (x-y)  under  nonnull  to  null  is  precisely  R.  However,  from  (2.3),  R  can  be 
expressed  as 


(2.4) 


R  *  E{q^  [  (c-<5)2+x’x^]/q^  [c2+xjx^] } 


where  E  stands  for  expectation  with  respect  to  c  having  the  pdf 
(2.5)  c  ~  q-j^ [c^x^ ]q2  [  (c-z;  2+y^y x  1  /  j  qi  [c2+x’x1  ]q2  [  (c-z)  2+y’yx ]dc 

•  — OO 


or 


:'s  -c.  l ' ■  ’  •.  •  -■  / 
Avjijob.*  i  U  y 

Viiat  j 


We  now  make  the  following  assumption 
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Assumption  2.1.  q,  ,q„e  Q  -  (q:  +  c,  for  every  a  >  0,  t>0}. 

1  2  q[c+t] 

It  then  follows  that  the  family  of  distributions  of  c  in  (2.5),  generated  by 
the  parameter  z,  keeping  x’x1  and  yjy  fixed,  has  MLR  property  in  c  by  the  assump¬ 
tion  q^  €  Q  and  hence  R  is  nondecreasing  in  z  by  the  assumption  q^  e  Q  (Lehmann 
(1959),  p.  74).  We  have,  therefore,  proved  the  following. 

Theorem  2.1.  For  testing  HqI  =  vs  Hi:  ®l  >  9  2  in  C^e  model  (2.2),  the 
test  with  a  critical  region  z  >  k  is  UMPI  in  the  class  of  conditional  level-a 
tests,  conditionally  given  the  two  ancillary  statistics  x^x^  and  yjy^*  whenever 
ql»q2  e  Q* 

Remark  2.1.  It  is  clear  from  (2.3)  that  when  both  q^  and  q^  are  normal, 
xjx^  and  y^y^  can  be  ignored  and  the  above  test  is  UMPI  unconditionally.  For 
nonnormal  q's,  sufficiency- invariance  reduces  (x,y,  x'x  ,  y{y..)  to  (x-y,  x^x.|,y!y1) 
and  the  above  conditional  argument  is  necessary  unless  n^n^-O.  This  means 
the  test  is  UMPI  when  there  is  only  one  observation  from  each  population  and 

q1»q2 c  Q* 

Remark  2.2.  If  qeQ  is  differentiable,  Q  is  equivalent  to 

(2.6)  Q  =  (q:  (-^,^2+-t))  +C,  V  t>0}. 

q(.c  +t) 

~  _yr  ^ 

An  example  of  a  q  e Q  is  provided  by  q(u)  *  e  ,  r  >  y.  It  may  be  remarked  that 
a  t-density  i  Q. 

Remark  2.3.  Our  attempt  to  derive  a  UMPI  or  even  a  LBI  test  when  q^  and  q^  involve 

an  unknown  common  scale  parameter  is  not  successful  unless  q^  and  q^  are  normal. 

It  is  possible  that  the  familiar  Fisher*s  t-test  for  this  problem  based  on 
1/2 

(x-y) /(xjx^+yjy^)  is  optimum  only  for  normal  densities. 
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Remark  2.4.  The  null  pdf  of  Z  is  given  by 


(2.7) 


<vd.  *  r  vc2H2«~>2>d' 

0  1  —CO 


where  q^(c  ) 


ni 

R  1 


q1(c2+x]!x1)dx1,  q2C(c-z)2)  = 


.n 


q2((c-z)  +yjy1)dy1.  Thus 


the  null  robustness  of  the  test  does  not  hold. 


3.  Test  of  Equality  of  Two  Scale  Parameters  (Non-Normal  Case) . 

Let  x^  be  * 1  random  vector  with  pdf  (with  respect  to  Lebesgue  measure) 

o.D  V*i>"ai  iqi((^i " yi!i} ' (?i * ui!i)/ai)  (i“1»2>» 


n. 


i  i  o 

where  e  R,  e^  *  (1,...,1)'  e  R  and  x^R  and  x2  e  R  are  independent.  Here 

2  2  2  2 

we  consider  the  problem  of  testing  H:  •  o2  vs  H^:  >a^.  Of  course,  if  x^’s 

are  both  normal,  then  the  F-test  is  UMPI.  We  shall  show  it  is  also  UMPI  in  the 
above  situation  under  a  mild  condition  on  q^,q2.  implies  an  optimality  ro¬ 

bustness  of  the  F-test.  Without  essential  loss  of  generality,  assume 
Then  the  problem  remains  invariant  under  the  transformation  x .  ■+  cx  (i-1,2), 
where  ceG  ■  (O,00).  Then  the  ratio  R  of  the  pdfs  of  a  maximal  invariant  under 
nonnull  to  null  is  given  by  R  ■  N/D,  where 


(3.2) 


N2N(o1,o2)  - 


r„2  "ni  2  ,  .2.,  n-1, 

,tni*lCTi  x;xi/oi)]c  dc 


,1/2, 


n*n^+n2  and  D  *  N(o^,a^).  In  (3.2),  transforming  c  into  c(x'x2)  / a 2»  N 

becomes 


(3.3) 


-n/2  n  .  fco 

NY  "  C?2?2}  Y  j  q1C-Ytc2)q2(c2)cn’3 


dc 


2,2 


where  y  -  °2^ai  and  t  *  xjx^/x2x2’  Note  D  *  N]_* 


Further  let  be  the  class  of 


functions  from  [0,®)  into  [0,®)  such  that  for  any  qe  0  ,  I  q(x'x)dx  ■  1  and 

”  pn 

q(gu)/q(u)  is  increasing  in  u  for  any  0  <6  <1.  Then  we  obtain 


2  2 

Theorem  3.1.  The  test  based  on  t  »  x’x./x’x.  >  k  is  UMPI  for  testing  a  *  cr.  vs 
2  2 

under  the  pdf’s  in  (3.1)  for  any  q^ e ,  q2  e  .  But  the  null  distribution 
of  t  depends  on  q^  and  q2> 

Proof.  Clearly  t  is  a  maximal  invariant  under  G.  Under  the  null  hypothesis,  the 
pdf  of  t  is  directly  shown  to  be 


(3.4) 


dP<  -  at'1"2*'2 
11 


2.  .  2X  n-1, 

q1(tc  )q2(c  )c  dc 


where  a  is  a  constant  depending  on  n.  and  n0  only  and  P*i  is  the  distribution  of 

1  4  H 

t  under  Hn:  a 2  ■  o2.  Hence  from  R  ■  dP^  /dP*  *  N  /N.  with  N  in  (3.3)  and 
0  12  ^  Hq  y  1  y 


(3.4),  the  pdf  of  t  under  (a^,o2)  given  by 

n,/2  i 


(3.5)  h(t |y)dt  -  RdP*  -  at(n-2)/2  y  1  [  q1(Ytc2)q2(c2)ctl”1dc]dt, 


so  h(t | y)/h(t [ 1)  ■  R.  Further  R  in  this  sense  is  expressed  as 

(3.6) 


2  2 

R  -  Y  1  E{q1(ycZ)/q1(cZ)] 


where  the  expectation  is  taken  with  respect  to  the  pdf 


(3.7) 


kt(.c2)  -  q1(c2)q2(.c2/t)cn"1/j00q1(c2)q2(c2/t)cn“1dc. 


Here  we  show  that  kt(c  )  with  t  >0  has  a  monotone  likelihood  ratio  property 

2.  ‘  ,  .2,  .  ,  /  2 , 


2  2 

in  c  .  In  fact,  since  q^  €  ,  k^  (c  )/k^(c2)  «q2(c2/t2)/q2(c2/t.^)  is  increasing  in  c 

2  22  2 

for  > ti*  The*efore,  since  q^(yc  )/q^(c  )  in  the  inside  of  the  expectation  of 


(3.6)  is  increasing  in  c  by  assumption,  its  expectation  is  increasing  in  t 


(see  Lehmann  (1959)  p.  74).  This  implies  R  in  (3.6)  is  an  increasing  function 
of  t.  Consequently  by  the  Neyman-Pearson  Lemma,  the  test  based  on  t  >  k  is  UMP. 
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Since  t  does  not  depend  on  q^’s,  it  is  UMPI  for  all  e  ,  q^  e  .  But  the 

null  distribution  of  t  in  (3.4)  cannot  be  free  from  q^'s.  This  proves  our  result. 

Remark  3.1.  If  “  0  is  not  assumed,  the  increasingness  of  q(Bu)/q(u)  is 

_  -  P  2 

simply  replaced  by  that  of  q(8u)/q(u),  where  q(u)  =  q(u+x  )dx. 

—00 

Remark  3.2.  A  crucial  distinction  of  our  result  from  Corollary  4.3  in  Kariya  (1981) 
lies  in  assuming  the  independence  of  x^  and  x^. 

Remark  3.3.  The  null  distribution  of  t  is  given  by  (3.5)  with  y  =  1.  Hence  for 
specific  q^  and  q^,  significance  points  of  t  are  obtained  from  it. 

Remark  3.4.  If  q  is  differentiable,  a  necessary  and  sufficient  condition  for 
q(Bu)/q(u)  to  be  increasing  in  u  for  0<B  <1  is  given  by  q'(u)  <0  and 

+  u.  This  condition  is  satisfied  for  a  multivariate  t-distribution  and  also 
for  more  general  normal  mixtures. 


4.  Test  of  Equality  of  Two  Location  Parameters  (Non-Exponential  Case). 


Let  X:  n  *  1  and  Y:  n  x 1  be  independent  random  vectors  with  pdfs  (with 

<«*  w 

respect  to  Lebestue  measure) 


>.l)  f0  (x)  -  qx[  l  (aci-ei)  ] ,  4  (y)  *  q2t  l  Cy^)  ] 

1-1  *  1 


x±  >_  6;l  ,  7i  ^  »X€^  ,  y  e  R  *  91  ’ 92  e 


The  problem  of  testing  H^:  9^  *  9^  vs  H^:  rema*ns  invariant  under  the 

group  G  of  translations  acting  as  x  -*x  +c,  y.  -*-y.+c,  i*l,...,n,  0  -*0  +c, 

till  11 

92^92+C’  -a’<c<"«  Let  xa)  *  min  x^y^  -  min  y±t  ^  *  I(xi“x( i))» 

n2  1 

t2  *  \  Applying  Wijsman's  (1967)  theorem,  the  ratio  R  of  the  pdfs 

of  a  maximal  invariant  T  *  (xqj-Yqj  ,  tl,t2^  unc*er  nonnull  to  null  is  obtained 


n, 


T  f  q-, £  I  (x.+c-e  )]q2[  l  (y,+c-e  )]dc 

'c>max(0  -x  ,0  -y  )  1  1 

<4'2>  ^ ^ - 

ti0  qx  t  I  (x^-Oj)  ]q2  [  ^(yi+c-ei)  ]dc 

Jc>max(ei-x(.1j,ei-y^1^)  1  1 


le^,(9l-^are2-ya^ 


qiItl+nl  (x(l)+C"0l)3q2ft2+n2(y(l)+C“02^dC 


•  c>oax[9  -x  9-T  ,VVVxCl)4^l)],2tt24V^l)4*'VJdc 
c_max  o x  x(1),01  y(1)j 


q, [t.+cjq-It  +  (c-z+6)  ]dc 

c>max(0,z-6)  “l 


q  [t  +c]q  [t  +  -=■  (c-z)  ]dc 
c>max(0,z)  nl 


q1[t1+c-6]q  [t  +  -=■  (c-z)  ]dc 
JC>max(6,z)  11  2  2  \ 

[  q.  [t  +c]q  [t  + -2-  (c-z)  ]dc 

' c>max(0,z)  11  1  1  nl 


•  ,“x(i)'7(i)-  {*9re 


To  derive  an  optimum  invariant  test  based  on  R,  note  that,  conditionally  given 
the  two  ancillary  statistics  t^  and  t^,  the  ratio  of  the  conditionalpdfs  of  z 
under  nonnull  to  null  is  precisely  R.  Moreover,  from  (4.2),  R  can  be  expressed 


(4.3) 


R  -  E{q1[(t1+c-0)}/q1[t1+c]}.I[c>(5] 


where  I  is  the  indicator  function  and  E  above  stands  for  expectation  with 
respect  to  c  having  the  pdf, 

”2  f  n2 

(4.4)  c  ~q  .[t  +c]q  +— -  (c-z)]/  q  [t  +c]q  [t  +—  (c-z)]dc,  c>max(0,z), 

1  *  c>max(0,z)  11  1  1  ni 
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We  now  make  the  following  assumption. 

Assumption  4.1.  ^i»q2  e  ^  =  ^q:  9(t+c-a)/q(t+c)  +  c  for  c  >  a,  Vt^O). 

It  then  follows  that  the  family  of  distributions  of  c  in  (4.4),  generated  by  z 
for  fixed  t^,t2>  bas  MLR  property  in  c  by  the  assumption  q2  €  <2.  and  hence  R  is 
nondecreasing  in  z  by  the  assumption  e Q,  [vide  Lehmann  (1959),  p.  74J.  This 
proves  the  following  result. 

Theorem  4.1.  For  testing  H^:  9^ = vs  H^:  9^ > 8^  in  the  model  (4.1),  the  test 
with  the  critical  region  z  >  k  is  UMPI  in  the  class  of  conditional  level-a  tests 
conditionally  given  the  two  ancillary  statistics  t^  and  t^,  whenever  q^,q2  e  2.* 
Remark  4.1.  It  is  clear  from  (4.2)  that  when  both  and  q^  are  exponential, 
t^  and  t^  can  be  ignored  and  the  above  test  is  UMPI  unconditionally.  For  non¬ 
exponential  q's,  sufficiency- invariance  reduces  (x^  ,y(i)  t0 

(x^^-y  jt^jt^)  and  the  above  conditional  argument  is  necessary  unless 
in  which  case  t^  and  t^  become  vacuous. 

Remark  4.2.  If  q  e  Q.  is  differentiable,  Q.  is  equivalent  to 

(4.5)  Q.  ■  (q;  cqf  Ct+c) /q(.t+c)  i  c,  Vt^O}. 

_  r 

An  example  of  a  q  e  Q.  is  provided  by  q(u)  =  e  ,  r  >  0. 

Remark  4.3.  The  null  pdf  of  z  is  given  by 

z  f  ~ 

(4.6)  dP  /dz  »  q  [c]q  [c-z]dc 

u  c>max(.0,z)  n2~l 

-  f00  nl  *  >  r<*>  2 

where  q^c)  =  q^+t^^  2  dt^  q2(c-z)  =  J  q2[c-z+t2]t2  dt2<  Thus  the 

null  robustness  of  the  test  does  not  hold. 


10 


5.  Test  of  Equality  of  Two  Scale  Parameters  (Non-Exponential  Case). 


Let  X:  n^ x 1  and  Y:  n^ *  1  be  independent  random  vectors  with  pdfs  (with 
respect  to  Lebesgue  measure) 


“•n  ‘  °i  V  l  <xr“i)/oi]’  fv2,o2(yJ  ‘  ff2  V  \(v i-u2>/021 

^  y  >  Vi  j  — .  oo  <  ^  0  • 

We  consider  the  problem  of  testing  H^:  *  a ^  vs  when  y^  and  are 

unknown.  The  problem  remains  invariant  under  the  group  G  of  transformations: 

Xi"*’aXi+Cl*  yi^ayi+c2*  ui'*‘avi+c1»  ^2  a^2+C2  *  cr2’”aa2’  ~“<c1»c2  <CD»  a>0* 

ni  n2 

Let  x(1)  »  min  x^  -  min  y±,  tfl(xi~x^)  ,  t2  -  I(yi"y(i))*  Applying 

Wijsman’s  (1967)  theorem,  the  ratio  of  the  pdfs  of  a  maximal  invariant  Z^t^/t^ 
under  nonnull  to  null  is  obtained  as 


(5.2)  R  = 


-n  f  at  +n  <ax  +c  -u  )  at  +n  (ay  )+c  -u  )  n  +„  -1 

dpZ  °2  j  ’i1 - q - Ji2[ - r - la  dadcldc 

H1  cl>“l~‘xfl)  ,1:2>“2~ayfn  ,a>0 _ _ 

dPZ.  ’  'n2  f  ,«  -ta  <a*(  ,+c  -a  •at24”2<Wm+c2’1,2),  "lV1.  .  . 

H0  Cj  |qx[ - 5^ - Jq2i - - 1«  dadcldc2 

cl>ul'alt(l)>c2>“2-ay(l)’a>0 


at  at  -(n,-l ) 


r.  at,  ,  ac2 

r.  ati  -  ato  -<n--l) 

joqi£~Jq2£— ]ct1  da 


,  q  (x  )  =  q  [  x+u  ]  du 

^0 


[V-.1V.U.  -  (a  -1 ) 

— -  .  6  * 

r  9, [az]q„[a]da 


.  „  02  ll 

,  6  *  —  ,  z  =  — 

al  C2 


q,  U6]q  [az  ]da  ,  > 

Jo  . 2  ’  ' 

"  "  ‘  *  '  —  ’  ■  "  0 

qi[a3q2£az_:L]<ia 
- (n_-l ) 

*  6  E[q1[a6]/q1Ia]] 

where  E  above  stands  for  expectation  with  respect  to  a  having  the  pdf 
(5.3)  a~q1[a]q2Iaz“'1]/J  q1[a]q2(az~1]da,  a  >  0. 

We  now  make  the  following  assumption. 

•m 

Assumption  5.1.  q^,q2  e  Q.  *  (qs  q[ad]/q[a]+  a  for  0  <  d  <  1}. 

It  then  follows  (vide  Lehmann  (1959),  p.  74)  that  the  family  (5.3)  has  MLR 
property  in  a  and  hence  R  is  nondecreasing  in  z.  Thus  we  have  proved  the 
following. 

Theorem  5.1.  For  testing  H^:  *  a 2  vs.  in  the  model  (5.1),  the  test  with 

the  critical  region  z  >  k  is  UMPI  whenever  q^,q2  e  2» 

Remark  5.1.  If  either  or  u2  or  both  are  known,  t^  and  t2  are  suitably  tfede- 

ti  "Hi 

fined  and  the  assumption  5.1  is  modified  accordingly.  1  2 

-(n2+3)/2  2 

Remark  5.2.  The  null  pdf  of  z  is  given  by  z  q  [a]q9[az  ]  a  da. 

J0  1  1 

Hence  for  specific  q^  and  q2  signigicance  points  of  z  are  obtained  from  it. 

Remark  5.3.  If  q  admits  a  derivative,  Q.  is  equivalent  to 

2.  *  (q:  aq'  (,a)/q(a)  +  a  >  0}. 

_  r 

M  -*U 

An  example  of  a  q  e Q  is  provided  by  q(u)  *  e  *  r > 0. 
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